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Abstract 

We systematically evaluate observables for hard exclusive electroproduction of 
real photons and compare them to experiment using a set of Generalized Parton 
Distributions (GPDs) whose parameters are constrained by Deeply Virtual Meson 
Production data, nucleon form factors and parton distributions. The Deeply Virtual 
Compton Scattering amplitudes are calculated to leading-twist accuracy and lead- 
ing order in QCD perturbation theory while the leptonic tensor is treated exactly, 
without any approximation. This study constitutes a check of the universality of 
the GPDs. We summarize all relevant details on the parametrizations of the GPDs 
and describe its use in the handbag approach of the aforementioned hard scatter- 
ing processes. We observe a good agreement between predictions and measurements 
of deeply virtual Compton scattering on a wide kinematic range, including most 
data from HI, ZEUS, HERMES, Hall A and CLAS collaborations for unpolarized 
and polarized targets when available. We also give predictions relevant for future 
experiments at COMPASS and JLab after the 12 GeV upgrade. 
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1 Introduction 



For the last 15 years, the handbag approach to hard exclusive leptoproduction of pho- 
tons (DVCS) 4 and mesons (DVMP) off protons has been extensively investigated both 
theoretically and experimentally. The handbag approach bases on factorization into hard 
(short-distance) partonic subprocesses and soft (long-distance) hadronic matrix elements 
[1, 2, 3, 4]. The latter are parametrized in terms of Generalized Parton Distributions 
(GPDs) [2, 5, 6]. The GPDs encode information on the longitudinal momentum distri- 
butions of the partons inside the proton as well as on the transverse localization of the 
partons [7, 8, 9]. The forward limits of some of the GPDs are the usual parton distributions 
(PDFs) and the lowest moments of the quark GPDs are related to the form factors of the 
proton. The GPDs give access to the total angular momenta of the partons making up the 
proton via Ji's sum rule [2]. Another important property of the GPDs is their universality, 
i.e. the same GPDs occur in DVCS as well as in DVMP although in different flavor com- 
binations. At present analytic methods to compute GPDs from QCD are lacking. Only 
lattice QCD provides numerical results on the lowest few moments of u and d quark GPDs 
[10, 11] for unphysical pion masses. There are also a number of models for GPDs available 
(e.g. Refs. [12, 13, 14, 15]) which can be confronted to experimental data on hard exclusive 
reactions in order to learn about GPDs. 

In the late nineties when the phenomenology of hard exclusive reactions commenced, 
estimates of observables for such reactions were made on the basis of simple ansaetze for the 
GPDs, see for instance Refs. [16]-[22]. As it turned out in the course of time, these GPD 
models were insufficient to account for the increasing amount of accurate data coming from 
HERMES and Jefferson Lab. More complex parametrizations of the GPDs were invented 
and utilized to analyze independently the data on DVCS [23, 24] as well as DVMP for light 
vector mesons [25, 26] and for pions [27, 28]. It should be noted that in DVCS frequently 
only reduced amplitudes, the so-called Compton Form Factors (CFF), have been extracted 
[29, 30] which represent valuable constraints on the GPDs. 

It is only recently that the universality property of GPDs has fully been exploited and 
a combined analysis of DVCS and DVMP carried through or a set of GPDs extracted from 
either DVCS or DVMP used to evaluate the other reaction. Thus, for instance, Meskaukas 
and Miiller [31] performed a combined analysis of the HERA data on DVCS and DVMP. 
These authors also used the GPD H advocated for in Refs. [25, 26] from an analysis of 
DVMP data, to compute DVCS observables for HERA kinematics. For a first attempt to 
compute DVCS along the same lines also for other kinematical regions see Miiller et al in 
Ref. [32]. The results look quite promising, no severe discrepancy has been observed by 
these authors. In the present article we are going to study systematically predictions for 
DVCS in an leading-order (LO), lowest-twist calculation of the DVCS amplitudes using 
GPDs determined in Refs. [25, 26, 28]. In contrast to earlier work, e.g. Ref. [21], the 
leptonic tensor is evaluated without any approximation as it is done in Ref. [33] implying 

4 Formally, Deeply Virtual Compton Scattering or DVCS refers only to the sub-process j*p — > -yp. 
However, DVCS is often used more loosely in the literature to name the photon leptoproduction process 
Ip — » Ipj used experimentally. 
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the inclusion of effects that are suppressed by powers of 1/Q in our analysis. The primary 
goal of our study is to examine how realistic the considered set of GPDs is in order to 
eventually improve it if necessary. 

The plan of the paper is as follows: in Sect. 2 we recapitulate the parametrization of the 
GPDs under scrutiny, and the theoretical description of the DVMP and DVCS processes 
in the handbag approach. In the second part, Sect. 3, we systematically compare model 
expectations to existing DVCS measurements. In Sect. 4 we give some predictions for 
observables which will be measured by the COMPASS, CLAS and Hall A collaborations 
in the near future. Finally, our summary and an outlook are presented in Sect. 5. 



2 Theoretical description 

2.1 The parametrization of the GPDs 

In this section we recapitulate the parametrization of the GPDs used in Refs. [25, 26, 28] 
and specify all the ingredients required to fit the DVMP data. The GPDs are functions 
of three variables, the usual invariant momentum transfer, t, the skewness defined as the 
ratio of light-cone plus components of the incoming (p) and outgoing {p') proton momenta 



' + 



(p + p 

In the generalized Bjorken regime of large Q 2 , large W but fixed x-q, it is related to 
Bjorken-x, Xb = Q 2 /(2p • q) by : 

where q is the momentum of the virtual photon and Q 2 its virtuality. The generalized 
Bjorken regime is defined by large Q 2 and large photon-proton cms energy W, but fixed 
Xb- The third variable, x, represents the average momentum fraction the emitted (k) and 
reabsorbed (k') partons carry with respect to the average proton momenta 

x - {k + k ' )+ (3) 

GPDs further depend on a factorization scale [5], which is usually taken as the photon 
virtuality unless specified otherwise. At LO, leading-twist accuracy, only the GPDs F = 
H, E, H, E contribute to DVCS which are characterized by the fact that the emitted and 
reabsorbed partons possess the same helicity. In Refs. [25, 26, 28] an integral representation 
of the GPDs is used 

fl rl-\p\ 

F\x, £,t)= dp d V 5(p + tr)-x) MP, V, t) + Di(x, t) Q(? - x 2 ) , (4) 
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where fa is a double distribution [5, 34] 5 and Di is the so-called D-term [35] which appears 
for the gluon and flavor-singlet quark combination of the GPDs H and E. The label % refers 
to specific quark flavors (or appropriate combinations) or to gluons. The D-terms only 
contribute to the real parts of the amplitudes. At small skewness however, the amplitude 
is dominated by the imaginary part. In Refs. [25, 26, 28] the D-terms are neglected. The 
advantage of the double-distribution ansatz for the GPDs is that polynomiality of the GPDs 
is automatically satisfied. For the GPDs H and E the D-term ensures the appearance of 
the highest power of the skewness in the Mellin moments of the GPDs. 

It is popular to write a double distribution as a product of a zero-skewness GPD and 
a weight function [34] 

w(p v) = r(2n * + 2) K 1 - \p\? -rfY 1 (5) 

2 2 ^+ 1 r 2 (n i + l) (l-|p|) 2ni+1 ' {) 
that generates the £ dependence of the GPD : 

fi(p,V,t) = F i (p^ = 0,t)w l (p,r ] ). (6) 

In Refs. [25, 26, 28] the parameter rii is taken as 1 for valence quarks and as 2 for sea 
quarks and gluons. The zero-skewness GPD is parametrized as the forward limit of that 
GPDs multiplied by an exponential in t 

F l (p,£ = 0,0 = F\p^ = 0,t = 0) exp (tp fi (p)) (7) 

The profile function, Pfi(p), is parametrized in a Regge-like manner 

Pfi(p) = -a' fi \np + b fi (8) 

where a' represents the slope of an appropriate Regge trajectory and b parametrizes the 
t dependence of its residue. In Ref. [36] a more complicated profile function for valence 
quarks has been proposed 

p fi (p) = {a' fi In l/p + b fi ) (1 - p) 3 + Afi p(l - pf (9) 

and exploited in an analysis of nucleon form factors. At small x the Regge-like part 
dominates while, for x — > 1, the term oc A^ takes the lead. It turned out that there is a 
correlation between x and t: The small (large) x behavior of the profile function controls 
the nucleon form factors at small (large) —t. For parametrizations like (7) without nodes 
except at the end-points, this correlation also holds for other moments of the GPDs and 
even for convolutions with hard scattering amplitudes. Thus, the Regge-like profile function 
is a sufficiently accurate approximation at small —t, the region we are interested in. It has 
also been shown in Refs. [37, 36] that, at large x the Regge-like profile function leads to 



5 The variables p and rj are usually denoted by j3 and a, respectively. However, we do not use this 
notation here in order to avoid a clash of notation. These latter symbols are already used for powers in 
the functional form. 
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an unphysically large distance between the struck parton and the cluster of the spectators. 
This distance provides an estimate of the size of the proton as a whole. 

The decomposition of the quark double distribution f q for the flavors q = u and d into 
/* al and f^ ca is done following the convention defined in Ref. [38] : 

/vai(P> V, t) = [f q (p, V, t) + e f f q (p, rj, t)] 6(p) , 

fL(p,v,t) = f q (p,V,t)e(-p)-e f f q (-p,r ] ,t)e(p). (10) 

where €f = +1 for F = H and E and —1 for H and E. For H and H and t — this 
prescription corresponds to the usual decomposition of parton distributions into valence 
quark and sea quark distributions, e.g. g va i = q — q. The GPDs respect the symmetry 
relations 

F°(-x,U) = e f F°(x,U), 

F^(-x,U) = -e f F s Ux,U), (11) 

and 

F^i-x^t) = 0, -l<x<-£. (12) 
For DVCS the C-parity even combination of the double distributions is required 

f^ + \p,V,t) = f q (p,rj : t)-e f f q (-p,rj : t). (13) 

The parameters appearing in the double distributions, in particular in Eq. (7) and 
Eq. (8), are fixed in an analysis of DVMP data in the kinematical region specified by 
^ < 0.1, Q 2 >3 GeV 2 , W>4 GeV and -t<0.6 GeV 2 . Cross section data are the only 
measurements available on a large Q 2 range (up to about 100 GeV 2 ) and they are dom- 
inated by contributions from the GPD H. Other GPDs significantly enter asymmetries 
which are only measured for 2<Q 2 <4 GeV 2 . Therefore evolution effects in this study 
are only sizeable for the GPD H and can safely be neglected for other GPDs [26, 28]. 
Since the scale dependence of H is available in [26], we still make use of it in our analysis 
of DVCS although it has practically no bearing on our results. It should be mentioned 
that in [25, 26] the evolution of the GPDs H and H is treated in an approximate way 
through the evolution of the PDFs in (7). A possible evolution of the profile function is 
ignored. At least for small skewness and small —t this approximation is reasonable as has 
been demonstrated in [26]. We follow this recipe which in any case is only of importance 
for the description of the DVCS cross section at HERA energies. It has been checked in 
Refs. [26, 28] that the valence quark GPDs are in agreement with the nucleon form factors 
at small —t and that all GPDs respect various positivity bounds [36, 39, 40]. At small —t 
there is also reasonable agreement between the moments of these GPDs and recent lattice 
results [10, 11]. But the GPD moments from lattice QCD have a flatter t dependence than 
those obtained from the GPDs we are discussing here and also flatter than nucleon form 
factor data exhibit. A possible explanation comes from the unphysical values of the pion 
mass used in present lattice simulations and the contamination from excited states [41]. 
Chiral extrapolations of the lattice moments have not yet been systematically performed, 
see for instance Ref. [42]. 
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2.1.1 Parametrization of the GPD H 



The GPD H is rather well determined since it controls the cross sections for electropro- 
duction of vector mesons for which a wealth of data is available. An advantage is that 
its forward limit occurring in Eq. (7), is a usual PDF. Therefore, only the parameters 
appearing in the profile functions have to be fixed. Linear Regge trajectories are assumed 

OLu = «/ii(0) + a' hi t with % = g, sea, val (14) 

As it is well-known the intercept controls the \ow-x behavior of the PDF [43]. A stan- 
dard Regge trajectory is assumed for the valence quarks, see Tab. 1. Since the sea- 
quark PDF is mainly driven by evolution for Q 2 > 4 GeV 2 it is furthermore assumed 
that ahsc&{t) = othg{t)- The gluon trajectory with an effective scale-dependent intercept 
ah g {0,Q 2 ), is directly seen in the HERA experiments [44, 45] (and references therein) and 
consequently fixed by these data. This trajectory is also quoted in Tab. 1. The trajecto- 
ries are accompanied by Regge residues assumed to have an exponential t dependence, see 
Eq. (8), with slopes taken as (m being the proton mass) 

fr/ival = , 

2 

b hg = b hsca = 2.58 GeV 2 + 0.25 GeV 2 In - . (15) 

It is convenient to expand the forward limits of H, the PDFs, in a power series of yfx 
(for p > 0): 

3 

H»(p,t = t = o) = p - 5 mi-p) 5 E c ^ /2 ' 

3=0 

3 

ir(j>,t = t = o) = (°)(i-p) 2 "* +1 5>,y/2 ? (i6) 

3=0 

where, with regard to the fact that the forward limit of H 9 is defined as pg(p), 

5 g = a hg (0) - 1 . (17) 

The expansion coefficients = Cij(Q 2 ) have been obtained from a fit to the CTEQ6M 
PDFs [46]; they are compiled in Tab. 1 too. 

The advantage of this expansion is twofold. First the integral (4) can be worked out 
term by term analytically leading to a corresponding expansion of the GPDs 

3 

Hi(x, £, t) = e b ^ CijHijix, £, t) . (18) 

3=0 

The integrals are given explicitly in Ref. [26]. Second, the term p~ Ss > in Eq. (16) 
guarantees that the longitudinal cross section on p° and electroproduction which is given 
by 

a L oc W 45 ° iQ ) (19) 
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gluon 


strange 


«val 




a(0) 


1.10 + 0.06 L-0.0027L 2 


ah n (0) 

rig v ) 


0.48 


0.48 


a' 


0.15 GeV~ 2 


0.15 GeV" 2 


0.9 GeV" 2 


0.9 GeV~ 2 


c 


2.23 + 0.362 L 


0.123 + 0.0003 L 


1.52 + 0.248 L 


0.76 + 0.248 L 


Cl 


5.43 - 7.00 L 


-0.327- 0.004 L 


2.88 - 0.940 L 


3.11 - 1.36 L 


c 2 


—34.0 + 22.5 L 


0.692 - 0.068 L 


—0.095 L 


—3.99 + 1.15 L 


c 3 


40.6 - 21.6 L 


—0.486 + 0.038 L 









Table 1: Parameters used for the GPD H, with L = In (Q 2 /Q 2 ) and Ql = 4 GeV 2 for the 
CTEQ6M PDF set. 



at fixed Q 2 and small xb (£), is in reasonable agreement with the HERA data [44, 45]. 
The data used in current PDF analyses, for instance Refs. [47, 48] or more recent ones 
[49, 50, 51], do not constrain the gluon PDF well for p<0.01. However, forcing the 
expansion of the gluon and sea PDFs to behave as p~ 5g at low p always leads to reasonable 
agreement of the DVMP cross section with the HERA experiments. With this prescription 
other sets of PDFs, e.g. Refs. [47, 48] , provide similar results as CTEQ6M. We stress that 
in all cases the expansions (16) are in good agreement with the original PDFs within their 
quoted errors. 

Finally, in accord with CTEQ6 analysis, the quark sea is simplified in Ref. [26] as 

TJU TTCl ITS 

-"sea — -"sea — ^s-" sea ' 

with k s = 1 + 0.68/(1 + 0.52 In Q 2 /Q 2 ), (20) 

where the Q 2 dependence of the flavor symmetry breaking factor k s was taken from the 
CTEQ6M PDFs. 



2.1.2 Parametrization of E 

Much less is known about E than for H. There is only the analysis of the Pauli form factor 
[36] which provides information on E for valence quarks. In addition there is an admittedly 
weak constraint from the asymmetries in electroproduction of p° mesons measured with a 
transversely polarized target [52, 53], for details see Sect. 2.2. 

The GPD E does not reduce to a PDF, the forward limit is not accessible in DIS. 
Therefore, the forward limit is to be fixed from exclusive experimental data as well. It is 
parametrized like the usual PDFs: 

^ai(P,e = 0,t = 0) = B-^l-^l + ^Kp—a-p)^, (21) 

where B(a, b) is Euler's beta function. The prefactor ensures the correct normalization of 
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the Pauli form factor at t — 0. Indeed the n — 1 moment 



-n0 



= f dpp n ~ 1 EUp,Z = t = 0) (22) 

reduces to K q which is the flavor-g contribution to the nucleon anomalous magnetic moment 
(k u = 1.67 , Kd — —2.03). The fits to the nucleon Pauli form factors performed in Ref. [36] 
fix the parameters specifying E for valence quarks: 

0^ = 4, = 5.6 . (23) 

In the spirit of the Regge model, the trajectory a eva i is taken to be the same as in H, see 
Tab. 1. The profile function (8) is evaluated with a' eval = a' hval and with slope parameters 
6 eva i taken to be zero. 

All this specifies the double distribution (6) for E^ al . The GPDs are then obtained 
from Eq. (4) where the factor (1 — p) 13 ^" 3 is expanded in a power series up to order 8 in 
order to perform the integration analytically. We stress that in the form factor analysis 
[36] the more complicated profile function (9) has been used while in the analysis of meson 
electroproduction [25, 26, 28] the small-p approximation (8) is adopted. There is an ongoing 
reanalysis of the form factor data [54]; preliminary results are close to those of Eq. (23). 
The uncertainties in the determination of E is reduced as compared to the results provided 
in Ref. [36]. 

The determination of E for gluons and sea quarks is still in its infancy. A rough 
estimate of these GPDs has been made in Ref. [55] along the lines proposed by Diehl and 
Kugler [56]. Again the double distribution construction is used and the forward limits of 
the gluonic and strange quark GPDs are parametrized as 

E s (p,£ = t = 0) = Nsp-^'^-pY*', 

E9(p,S = t = 0) = N g p- s «(l-p)te<>. (24) 

Of course the same Regge trajectory as for H is used, see Tab. 1 and Eq. (17). The lack of 
detailed information forces the assumption of a flavor-symmetric sea. The powers of the 
large p behavior are set to the following values (variant 3 of Tab. 1 in Ref. [55]): 

P s = 7, P 9 = 6. (25) 

The integer powers allow to solve the integral (4) analytically. The slope of the residues 
are taken as : 

b eg = b ea = 0.9b hg , (26) 

see Eq. (15). In Ref. [55] the normalization, N s , of E s is fixed from saturating a positivity 
bound [36] for a certain range of p. Since the bound is quadratic the sign of N s is not 
fixed. Therefore, we have to consider the two cases N s = ±0.155. 

The normalization of E 9 is fixed by using a sum rule for the second moments of E [56] 



20 ~~ c 20 e 20 



Con — 



2 E e 20 (27) 



i=u,d,s 



8 



The valence quark contribution to this sum rule is very small. Hence, the gluon and sea 
quark moments almost cancel each other. In any case, the sum rules allows to fix the 
normalization N g for given N s . 



2.1.3 The GPD H 

The forward limit of H reduces to the polarized PDF for which in Refs. [26, 28] the 
Bliimlein-Bottcher results [57] are taken. To fix the parameters of H only the HERMES 
data on the cross sections and the target asymmetries for n + electroproduction [58, 59] 
are at disposal. Therefore, H is determined only for the valence quarks, H SC3u and H 9 are 
neglected. 

Analogously to Eq. (16) if* al i s expanded in a power series 



HUP, e = * = 0) = ife A q p- a ~^ (i _ p f c q3 p> , (28) 

3=0 

where q = u,d and A q is a normalization factor. The GPDs are constrained by the lowest 
moments at £ = t — 

Vq = [ dpH^(p,t = t = 0) (29) 
Jo 

which are known from F and D values and /3-decay constants in flavor SU(3) 

r] u = 0.926 ±0.014, rj d = -0.341 ±0.018. (30) 
This normalization is guaranteed by the factor 

(2-^(1-^) 



A" 1 = B(l-a hq ,4) 



. . . 1 - a hq , - 
CqO + Cql Z r C q2 

5 - a~ hq 



(6-tt^)(5-a^) 



hq) 



(31) 



The expansion coefficients are compiled in Tab. 2 as well as the other parameters specifying 
this GPD. 



2.1.4 The GPD E 

Only exclusive ir + electroproduction data constrain this GPD in the DVMP analysis. 
Therefore, as for H, it can only be fixed for valence quarks. As is well-known it con- 
sists of two parts the pion-pole contribution and a non-pole one. The pole contribution 
reads [60, 61] 

EUe = -eU = e(M < ^ (^) , (32) 

where Fp is the pseudoscalar from factor of the nucleon. With the help of PCAC and the 
Goldberger-Treiman relation the pole contribution to the pseudoscalar form factor can be 
written as 

rp / 1\ f 2y/2gnNNF nNN (t) 

F P (t) = -m N f v ^ . (33) 

t — mi 
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dval 


a(0) 


0.48 


0.48 


a' 


0.45 GeV~ 2 


0.45 GeV" 2 


h 










0.170+0.03 L 


-0.320-0.040 L 


Cl 


1.340-0.02 L 


-1.427-0.176 L 


c 2 


0.120-0.40 L 


0.692-0.068 L 



Table 2: Parameters used for the GPD H. Evolution is parametrized through the variable 
L = m(Q 2 /Q 2 ) with Q 2 = 4 GeV 2 . 



Here m N (m n ) is the mass of the nucleoli (pion), g w NN (= 13.1) is the pion- nucleoli coupling 
constant, f n is the pion decay constant and is the pion's distribution amplitude taken 
as: 

^(r) = 6r(l - r) [1 + a 2 Cf{2r - 1)] . (34) 

The Gegenbauer coefficient is taken as a 2 = 0.22 at the initial scale Q 2 = 4 GeV 2 . This 
value for a 2 is conform with the sharp rise of the transition form factor with Q 2 as is 
observed by the Babar collaboration [62]. However, this behavior of the transition from 
factor is not seen by the BELLE collaboration [63]. The BELLE data are compatible with 
a pion distribution amplitude close to the asymptotic form 6r(l — r). Since E plays a 
minor role in DVCS we keep the value 0.22 for a 2 ; it practically has no bearing on our 
predictions for DVCS. 

Finally, F wNN is the form factor of the pion-nucleon vertex and is parametrized in 
Ref. [28] as 

A 2 ^ Ul^ 

FnNN = — 72 rf (35) 

with Atv = 0.44 GeV. 

The non-pole contribution for which there is only a weak evidence in the data, is 
parameterized in the same spirit as the other GPDs. The forward limit is parametrized as 
in Refs. [28, 64] 

EU P ,Z = t = 0) = Nlp^){l-pf. (36) 

Flavor independence of the Regge trajectory and the slope of its residue is assumed. The 
parameters for E are compiled in Tab. 3. 

2.2 Description of DVMP 

The particular variant of the handbag approach used to extract GPDs from meson elec- 
troproduction in the kinematical region of large Q 2 and large W but small xb( ^0.2) and 
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«e(0) 


<4 


h (GeV- 2 ) 


N? 


TVf 


0.48 


0.45 


0.9 


14.0 


4.0 



Table 3: Parameters used for the GPD E. Evolution is ignored. 



small invariant momentum transfer —t, is described in some detail in [25, 26]. Here, only 
the basic facts are reviewed. As an example, let us examine the helicity amplitudes for 
the asymptotically leading transitions from longitudinally polarized photons to likewise 
polarized p° mesons, 7£p — > p° L p : 

e i 

M -,o + = + e^ty - ejft} , (37) 

within the handbag approach. The helicities of the protons are labeled by their signs and 
e q denotes the charge of the quark with flavor q in units of the positron charge |e|. To the 
proton helicity-non-flip amplitude the GPDs contribute in the combination 

H cS = H--^-E. (38) 
1 - ? 

The terms Ty in Eq. (37) denote convolutions of subprocess amplitudes and GPDs F 
(=H,E): 

Hit, *, Q 2 ) = Y.j dx ^>w( x > Z> Q 2 > 1 = °) ^(x, £, t) (39) 

A Jx i 

where % = g,q, x g = and x q = —1. The subprocess amplitude A for partonic helicity 
A is to be calculated perturbatively. Since E is of the same order as H in absolute value 
(see the discussion in Sec. 2.1), one can approximate H e g by H in the small skewness 
region. Moreover, for small —t' the helicity-flip amplitude can also be neglected in the 
cross sections for vector mesons ( an exception is p + production) which is therefore only 
sensitive to H. 

As is well-known, in collinear factorization the handbag result for the integrated cross 
section scales as 1/Q 6 at fixed xb- However, for the kinematics accessible to current 
experiments, the data are in conflict with this prediction. This can be seen from Fig. 1 
- the recent HI data [44] on p° production only drop as ~ 1/Q 4 - The theoretical 1/Q 6 
behavior of the cross section is modified by logs of Q 2 generated by the evolution of the 
GPDs; they diminish the discrepancy between theory and experiment. Experiment also 
tells us that the transverse cross section, <tt, is not small; the ratio R = (Jl/pt, also shown 
in Fig. 1, is rather small for experimentally accessible values of Q 2 . Combining the data 
on the ratio with the unseparated cross section it becomes clear that the longitudinal cross 
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10 3 









10 2 : 


1 




t 




a, 

# 






10 1 : 







10" 




Q 2 [GeV 



Q 2 [GeV 2 ] 



Figure 1: Left: The cross section for p° electroproduction vs. Q 2 at xb — 0.002. Data are 
taken from Ref. [44] and compared to a power-law fit. Right: The ratio of <tl and gt for 
p° production vs. Q 2 at W = 90 GeV. Data taken from Refs. [44, 45]. The figure is taken 
from Ref. [26] where also further references to data can be found. 



section approximately falls off as ~ 1/Q 4 too. This fact implies a marked overestimate 
of the longitudinal cross section at Q 2 ~ 4 GeV 2 in the collinear approximation if it is 
evaluated from GPDs of the type discussed in Sect. 2.1.1. Agreement with experiment 
is however achieved in collinear approximation for an alternative parameterization of the 
GPD H proposed in Ref. [31]. The evolution of this GPD produces much larger logs of Q 2 
than the GPD described in this article. 

In view of this situation, power corrections are added to the LO subprocess amplitudes 
in Refs. [25, 26]. These power corrections are calculated within the modified perturbative 
approach [65] in which quark transverse degrees of freedom are kept and gluon radiation is 
taken into account. The latter has been calculated in the form of a Sudakov factor to next- 
to-leading-log approximation using resummation techniques and having recourse to the 
renormalization group [65]. For consistency, allowance is to be made for meson light-cone 
wave functions instead of distribution amplitudes. Due to these power corrections the con- 
volutions J-y also depend on Q 2 . The modified perturbative approach is designed in such 
a way that asymptotically the collinear result for the subprocess amplitudes emerges. It is 
to be stressed that, in contrast to the situation at the mesonic vertex, the partons entering 
the subprocess are treated as being emitted and reabsorbed by the proton collinearly. 

The above described approach can also be applied to vector mesons other than the p° 
and can be generalized to the asymptotically suppressed amplitudes for 7^ — > Vr tran- 
sitions 6 . In collinear approximation these amplitudes are infrared singular but in the 
modified perturbative approach in the propagators regularizes the singular integrals. 
It should be mentioned that the described approach bears similarities to the color dipole 

6 Another method to treat — > Vt transitions has been proposed in Ref. [66] 
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Figure 2: The longitudinal cross section of p° electroproduction vs. W at Q 2 = 4 GeV 2 . 
The handbag result is shown as a solid line; the shadowed band represents the uncertainties 
of this result. For references to the data see [26, 68]. Figure taken from [68]. 



model, see Ref. [67] and references therein. The described approach can also be applied to 
electroproduction of pseudoscalar mesons [28, 64] where one learns about the valence-quark 
components of GPDs H and E (see Sections 2.1.3 and 2.1.4) as well as on some of the 
transversity GPDs. However, the latter do not contribute to DVCS. 

The available data from HERMES, COMPASS, E665, HI and ZEUS on cross sections 
and spin density matrix elements for p° and <fi electroproduction have been analyzed in Refs. 
[25, 26]. The data cover a large range of kinematics: Q 2 varies between 3 and 100 GeV 2 and 
W between 5 and 180 GeV. In Fig. 1 the ratio o^I^t is shown for p° electroproduction 
in order to demonstrate that a fair description of the amplitude for 7^ — > Vr is also 
achieved. In Fig. 2 the longitudinal cross section for p° production is shown versus W 
at Q 2 = 4 GeV 2 . As the inspection of the figure reveals, a good description of all the 
1ow-xb data has been achieved for W > 4 GeV. For smaller W the handbag results deviate 
from experiment; at W = 2 GeV theory and experiment deviate by orders of magnitude. 
For cf) production, on the other hand, the handbag seems to work even at W ~ 2 GeV. 
More results and references to the experimental data can be found in Refs. [26, 68]. The 
HERMES cross section data on 7r + electroproduction data [58] as well as various single 
spin asymmetries are analyzed in Refs. [28, 64]. 

In the mentioned kinematical region, information on the GPD E can only be extracted 
from the target asymmetry 

r ^H w ] , (40) 



which is obtained through the sin {<p — 4>s) harmonic of the electroproduction asymmetry 
measured with a transversally polarized target. Here, <fi is the azimuthal angle between 
the lepton and the hadron planes and 4>s specifies the orientation of the target spin vector 
with respect to the lepton plane (in the Trento convention [69]). The convolutions Sy and 
TLy in Eq. (40) are defined in Eq. (39). It is shown in Ref. [55] that the GPD E described 
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in Sect. 2.1.2 in combination with the rather well-known GPD H, provides results in 
agreement with recent p° data on this asymmetry from the HERMES [52] and COMPASS 
[53] collaborations. It is to be stressed that for p° production the sin (0 — <f>s) harmonic 
of the cross section for a transversally polarized target essentially probes E for valence 
quarks; the gluon and sea-quark contributions cancel each other to a large extent due to 
the sum rule (27) which implies that the second moments of E 9 and E sea have about the 
same strength but opposite sign. For the simple parametrization (24) with no nodes except 
at the end-points, this property of the second moments transfers to other moments of E 
to a certain degree and in particular to the convolutions. In accord with this argument, 
Al^^~^ for electroproduction of the mesons is predicted to be about zero in agreement 
with preliminary HERMES data [70]. 

2.3 DVCS in the handbag approach 

In the following we use the aforementioned GPDs constrained from DVMP, nucleon form 
factors and partons distributions to evaluate Ip — > Ipr) observables and compare to mea- 
surements. As for DVMP we will analyze this process in the generalized Bjorken regime 
of large Q 2 , large W but fixed xb- 

Leptoproduction of photons is complicated since, besides the DVCS contribution, j*p — > 
7P, there is also the Bethe- Heftier (BH) contribution. The square of the Ip — > Ip^j amplitude 
Mip^^ therefore falls into three parts : 

\Mi p ^i m \ 2 = \Mbh\ 2 + Mi + |.Mdvcs| 2 - (41) 

These parts readily correspond to the squared amplitudes of the BH and DVCS processes 
and their interference. In the one-photon-exchange approximation of QED, the three terms 
|-Mbh| 2 , |A^dvcs| 2 an d Mi in (41) have the following harmonic structure in 0, the az- 
imuthal angle of the outgoing photon with regard to the leptonic plane (in the Trento 
convention [69]) : 

1 1 3 

I " Mbh|2 K mp7oo^^ KHc ° sM) + S " Hsin(n0)] ' 
IMwcsl 2 oc ^[c° VCS cos(n0) + S ° vcs sin(n0)] , 

?1=0 

1 1 3 

Ml K W\ P(coscb) ^ K cos ( n ^) + 4sfn(n0)] , (42) 

' ' ^> 71=0 

where P(cos0) comes from the BH lepton propagators. Although there are only harmonics 
up to the maximal order 3 in the sums, the additional cos <fi dependence from the lepton 
propagators generates in principle an infinite series of harmonics for the BH and inter- 
ference terms. A more detailed harmonic structure taking into account beam and target 
polarizations, can be found for instance in [71]. For transverse target polarization the 
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harmonic series also depends on the angle 4>s which specifies the orientation of the target 
spin vector. An harmonic analysis of Ip — > Ip^y allows for an examination of the GPDs 
[72]. Detailed analytic expressions describing this harmonic structure were published in 
Ref. [21]. They involve CFFs, which are integrals of GPDs over the momentum fraction x 
with a hard scattering kernel. At LO the CFFs read : 



dx 



e 2 u F u + e 2 d F d + e 2 F* 



1 1 



x — is £ + x — ie 



(43) 



where e/ is defined after Eq. (10). The CFFs, the analogues of the convolutions (39) for 
DVMP, are complex functions due to the singularity at x = ±£ in the integration domain. 

In Ref. [21] the electroproduction of photons was evaluated to leading and subleading 
order in an 1/Q expansion. Since a great wealth of data in the valence region [73, 74] involve 
not-so- large values of Q 2 , the impact of this approximation is not negligible when comparing 
theoretical expectations and measurements. In 2008 Guichon and Vanderhaeghen [75] 
elaborated on their previous numerical computations [60, 76, 77] of the ep — > ep'j process 
to establish analytic expressions of the cross section for all polarizations of the proton target 
and of the incoming lepton. In these formulas the leptonic tensor is treated exactly, i. e. the 
involved kinematic terms are kept with their full Q 2 dependence : they are not expanded 
as power in 1/Q. These analytic expressions were implemented into a ROOT/C++ code 
[30] and it was checked in Ref. [78] that they are completely equivalent to the expressions 
of cross sections for all polarizations of the proton target and of the incoming lepton, 
presented in Refs. [60, 76, 77]. Later Belitsky and Mueller extended their earlier work [21] 
by removing the approximations done in the 1/Q expansion of the leptonic tensor [79, 33]. 
Note that the case of the transverse target polarization is not treated in these references. 

In the present work, we use the theoretical framework of Guichon and Vanderhaeghen 
which provides a complete and accurate set of formulas that encompass all types of po- 
larizations, including the case of a transverse target. Our approach is consistent with the 
so-called BM formalism [33]. 

Further improvements of this approach have been considered: the next-to-leading order 
(NLO) kernels have been calculated in Refs. [80, 81, 82, 83, 84]. Soft-collinear resummation 
formulas have been derived recently [85]. Finite-t and target-mass corrections to DVCS 
have recently been investigated in the OPE framework in Refs. [86, 87] (and references 
therein). Twist-3 effects in DVCS have also been studied in Refs. [88, 89, 90]. The phe- 
nomenological implications of all these corrections for DVCS observables remain to be 
studied. This is beyond the scope of the present article. 



3 Comparison to DVCS data 

The Ip — > Ipj cross section on an unpolarized target for a given beam charge, e; in units 
of the positron charge and beam helicity hi/ 2 can be written as : 

da h " ei (<j)) = rfauu(0) [1 + M lu,dvcs(0) + e^i^LU,i(0) + e;A c (0)] , (44) 
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where only the dependence of the observables is shown. If both longitudinally polarized 
positively and negatively charged beams are available, the asymmetries in Eq. (44) can be 
isolated, as is the case for a large part of HERMES data. Thus, for instance the beam 
charge asymmetry is obtained from the combination : 



M<t>) 



1 



4<i<7uu(</>) 



(da^ + da"-) - (da^ + da^) 



(45) 



From analogous combinations, one obtains the two beam spin asymmetries Am,i and 
Alu,dvcs'- 



AluM 
^4lu,dvcs(0) 



4dcxuu(0) 
1 

Ada vv ((j)) 



(da^ - da*-) - (da^ - da*~ 
(da^ - da^) + (da^ - da^)j . 



(46) 
(47) 



If an experiment only has access to one value of ei such as in Jefferson Lab, the asym- 
metries defined in Eq. (44) cannot be isolated and one can only measure the beam spin 
asymmetry A^u which depends on the charge-spin cross section as follows : 



KM = 



da — da* 
da^ + da* 



(48) 



where we use the familiar notation of labeling the charge-spin cross section by the sign of 
the beam charge e\ and an arrow — > (■<— ) for the helicity plus (minus). One can check that 
A^u can be written as a function of the spin and charge asymmetries defined in Eq. (44) : 



^U(0) = 



a A 



LU,R 



A 



LU,DVCS 



(0) 



l + e z A c (0) 



(49) 



The case of longitudinally polarized target observables is simpler, due to the fact that 
there are no data with both varying longitudinal target polarization and beam charge. 
Therefore, experiments measured the target longitudinal spin asymmetry which reads : 



^l(0) = 



[da^ + da~ 



[da' 



+ da^ 



[da*~^ + da^^] + [da*~^ + da 



(50) 



where the double arrows -<= (=^) refer to the target polarization state parallel (anti-parallel) 
to the beam momentum. The double longitudinal target spin asymmetry is defined in a 
similar fashion : 

e i e i e i e i 

[da^^ + da*~^\ — [da*~^ + da^^] 



A 



ll(0) = 



[da^ l=> + da* 



[da*~ l=> + da^ 



(51) 
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The HERMES collaboration also had access to a transversally polarized target with 
both electrons and positrons. They therefore were able to measure two types of observ- 
ables : 



-4ut,i(0, 0s) = 

da + ((f), 0s) — da + ((j) : 05 + 7r) + da~((f), 0s) — da~((f), 0s + 7r) 
da+ (0, S ) - da + (0, 0s + n) + da~ (0, <f> 3 ) - da~ (0, 0s + n) 

^UT,DVCs(0, 0s) = 

da + (0, 0g) - da + (0, 0g + tt) - dcr" (0, 5 ) + (icr~ (0, 0g + tt) 
da+ (0, 0g) - da+ (0, S + 7r) + dcr" (0, 0s) - do- (0, 5 + tt) 



,(52) 



.(53) 



Finally, it is worth mentioning that the HERMES collaboration usually does not publish 
the 0-dependence of the asymmetries but rather chose to extract harmonics out of their 
asymmetries. For instance, in the case of the beam charge asymmetry Ac, the cos(n0) 
harmonics are extracted using the following formula : 



/•27T 

^cosM) = N / d( p Ac ^ cos ( n 0) ? (54) 
Jo 

where the normalization factor iV is 1/2tx in the case n = and 1/n for n > 1. Since 
HERMES has not measured numerators (V) and denominators (S) of the asymmetries in 
Eqs. (45) - (53) separately the projection (54) is the best approximation to the Fourier 
coefficients 

1 j^d(f)V{(f)) cos(n0) 

2 | o 27r #<S(0) 

HERMES can do. The latter Fourier coefficients are closest related to the CFFs. 

The coefficients of the various Fourier harmonics occurring in Eq. (42) provide infor- 
mation about CFFs, or equivalently GPDs, either in the interference with the real Bethe- 
Heitler amplitude or from the DVCS process. Through the measurement of cross sections 
or the various asymmetries described in Eqs. (45)-(53), one can put constraints on different 
combination of GPDs. As an illustration we quote the leading-twist, LO pQCD connection 
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between observables and CFF for some asymmetries [21, 33, 71] 



A c ° s<t> oc Re 



F 1 H + £(F 1 + F, 
~F X U + i{F x + F 2 )U - J- 



4m 2 



F 2 S 



A ul,i K Im 

oc Re 

,cos(0<« n 

^LL,DVCS ^ ^ 



jsin(<^.-0 s ) 
^■UT,DVCS 

^sin cos* ^ Im 



~m + F 2 )( W + ) + Fl H - aj^F, + 4^F 2 )E_ 

4(1 - ^ 2 ) + me) - 4^ 2 + £%* + + en*) 

Im(W£*) -{Im(Wr) 



.^(F,M-F 1 £)+?(F 1 + JLf 1 )(« + £) 



(56) 



In the following subsections, we compare predictions on DVCS evaluated from the GPDs 
described in Sect. 2.1 to the available data in different kinematic domains : from small 
x B (HI, ZEUS) [91]-[94] to intermediate x B (HERMES) [95]-[98] and large x B (Jefferson 
Lab Hall A and CLAS) [73, 74, 99]. The HERA HI and ZEUS collaborations measured 
the DVCS (7*p — > 7p) cross section at large W and Q 2 but small x B , which are mostly- 
sensitive to the imaginary part of H. HERMES measured and published a large number 
of asymmetries with different beam charge, helicity states and target polarization. These 
data lie in the range 0.05 < x B < 0.25 and are sensitive to the real and imaginary parts 
of the CFFs "H, S and H. The data sets from Jefferson Lab either cover a wide kinematic 
range (0.11 < x B < 0.58) [74] or are highly precise on a restricted kinematic domain [73]. 
In the first case beam spin asymmetries are measured, and in the second helicity-dependent 
and independent cross sections. The dependence of all these observables to the CFFs are 
compiled in Tab. 4 where for each experiment, we have chosen the typical kinematics listed 
in Tab. 5. The BH amplitude is evaluated from the Kelly parametrization of nucleon form 
factors [100]. We would like to stress once again that no DVCS data has been used in order 
to fix the GPD parameters. Observables which vanish to the accuracy we are working and 
for which the data are compatible with zero within errors will not be discussed in the 
following. The propagation of PDF errors to the GPDs H and H is the only source of 
uncertainties we have considered so far. They will be shown as shadowed bands around 
the predictions on all the following figures. 
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Experiment 


Observable 


Normalized CFF dependence 




a cos 0<p 


Ren + 0.06Re£ + 0.24Re% 




A COS <fr 


ReV. + 0.05Re£ + 0.15Re^ 




a sin <f> 


Imn + 0.05Im£ + 0.12Im^ 




^+,sin0 


lm.n + O.lOIm^ + 0.01Im£ 


HERMES 


a 1 G1V1 V/Tl 


Imn - 0.97Im% + 0.49Im£ - 0.03Im£ 




i —1— rriQ C\A\ 

A 1^ i tUB ijljiy 


1 + 0.05Re% + O.OlReT/ 




^ -1- COS (7) 


1 + 0.79Re"H + O.llIm'H 




,s'm((f)-<j>s) 
^UT.DVCS 


T 1 / V) O T C T ) 1 1 

ImnRec- — LmtReH 




a sin(<j)—(ps) cos <j> 


lm/i — U.oolmc — U.lzlm/i 




a — .sin <b 


IrnH + 0.06Im£ + 0.21Im% 


L>LAb 


< -,sin</J 


lm/i + U.lzlm/i + U.U41mc 




x — ,sin 20 


Im% - 0.79Im% + 0.30Im£ - 0.05Im£ 




A(T sin0 


Inift + 0.07Im£ + 0.47Im% 


HALL A 


£j-cos 00 


1 + 0.05Re% + o.oo7nn* 




^-cos <j!> 


1 + 0.12Re% + 0.05Re% 


HERA 


0"DVCS 


nn* + omss* + nn* 



Table 4: Dependence of the observables on the CFF at the kinematics specified in Tab. 5. The 
coefficients in front of the CFF are normalized to the largest one, but only relative coefficients 
larger than 1% are kept, except for the Hall A cross section, where we also show the term quadratic 
to the CFF T~L since it contributes significantly. In order to get simple CFF dependences for this 
table, the unpolarized cross section in the denominator of the asymmetries is approximated by 
the term c^ H /P(cos <p) in Eq. (42). 



Experiment 


Kinematics 


x B 


Q 2 [GeV 2 ] 


t [GeV 2 ] 


HERMES 


0.09 


2.50 


-0.12 


CLAS 


0.19 


1.25 


-0.19 


HALL A 


0.36 


2.30 


-0.23 


HERA 


0.001 


8.00 


-0.30 



Table 5: Typical kinematics used in Tab. 4 for various experiments. 
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Figure 3: The cos 00 and cos <p harmonics of the beam charge asymmetry at the kinematical 
setting x B ~ 0.097 and Q 2 ~ 2.51 GeV 2 . Data are taken from HERMES [97] -Tab. 6. Our 
results, shown as solid lines with the shaded areas as the error bands, are evaluated at the 
kinematics specified in Ref. [97] -Tab. 6, and joined by straight lines to guide the eyes. 



3.1 Beam charge asymmetry 

We begin the discussion with the beam charge asymmetry, Ac, as measured by the HER- 
MES collaboration during 1996-2007, a result which has been recently updated in Ref. [97]. 
This observable is generated by the BH-DVCS interference term (see Eq. (56)) which, to 
leading-twist accuracy, feeds the cosn^ harmonics (n = 0,1). As inspection of Tab. 4 
reveals, both the harmonics depend mostly on the real part of the CFF TL. In Fig. 3, our 
predictions are compared to the HERMES data [97]. Very good agreement can be seen for 
both the harmonics, demonstrating that in this kinematical range, the real part of % has 
the right magnitude. 

3.2 Beam spin asymmetry 

As can be seen from Eq. (56) the sin0 harmonic of the beam spin asymmetry, Amj, 
depends on the same combination of electromagnetic form factors and CFFs as the beam 
charge asymmetry As we already mentioned this combination is dominated by the GPD 
H, see Tab. 4. Since the electromagnetic form factors are real, the imaginary part of the 
CFF 7i is required which, to leading-order of perturbative QCD (see (43)), is given by the 
GPD at the cross-over line x = £, i.e. essentially probes the combination 



t [GeV 2 ] 



elH u (C, & t) + e 2 d H d {t £, t) + e 2 s H s (£, £, t) . 



(57) 
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Our results for the beam spin asymmetry A s ^fj are shown in Fig. 4-left and compared to 
the HERMES data [97]. The agreement between predictions and data is not as good in this 
case, our results differ by about 40% (~ 0.1 in absolute value) from experiment. Recently 
the HERMES collaboration has published data on the sin0 harmonic of the beam spin 
asymmetry using a recoil detector and a positron beam [98]. In this experiment all three 
final state particles are detected and therefore the resonant background severely reduced. 
In so far the recoil data are closer to the exclusive process Ip — > Ipj to which our theory 
applies. The data were taken at about the same average values of xb and Q 2 . In order 
to compare to recoil data, we computed A\ v using Eq. (49) with Alu,i an d Ac from the 
non-recoil data and Alu,dvcs = (exact at twist 2 and in agreement with experimental 
results from Ref. [97]). Then the sin0 coefficient is : 



.4 



+ sm<; 
LU 



a sin cf> 
n LU,I 



l + A 



(58) 



c 



On the right hand side of Fig. 4 we therefore show both from the non-recoil and 

the recoil data. We observe that the recoil data are significantly larger in absolute value, 
yielding very good agreement with our predictions. Similar effects for other DVCS ob- 
servables may occur but with the exception of the beam spin asymmetry, there are no 
measurements with the recoil detector available. The effect of the resonant background in 
other observables is unknown. Note that A^^ vanishes for forward scattering, t = t m[n . 
The trend towards zero is however only visible for t of order t min = — 4m 2 £ 2 / (1 — £ 2 ) which 
is very small, about -0.02 GeV 2 for HERMES kinematics. 




t [GeV ] 



Figure 4: Left plot: ^4lui as a function of —t measured by the HERMES collaboration 
[97]-Tab. 5. Right plot: A^ m( ^ versus —t obtained from the non-recoil data on ^luI an< ^ 
^coso</> measurec [ by the HERMES collaboration [97] -Tab. 5 and 6 (solid circles, see text 
for details) and the more recent recoil data [98] (open squares). For other notations and 
the values of the averaged kinematic variables, refer to Fig. 3. 



The CLAS collaboration published accurate data on the beam helicity asymmetry in a 
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large kinematical range [74]. The data analysis required the detection of the full (e,p, 7) 
final state, much like the HERMES recoil data. A selected sample of the 0-dependent 
asymmetry A^ v (<f)) is shown on Fig. 5. While the dominance of the sin0 term is clearly 
visible in the data one notices a discrepancy between the data and our prediction of the 
same order as for the non-recoil HERMES beam spin asymmetry, see above. We repeat - 
our GPDs are optimized for small skewness (see discussion in Sect. 2.1 and Sect. 2.2) and 
one therefore cannot expect perfect agreement with large-skewness experiment. 



^LU(^) 





-fx 


I \: 






(ir B )=0.13 V_/ 
(Q 2 ) = 1.17 GeV 2 


(x B ) = 0.18 M { / 
(Q 2 ) = 1.37 GeV 2 \-S 


' \ 

(x B ) = 0.25 \ I / 
(Q 2 ) = 1.69 GeV 2 \J/ 


- XI 1/ 

<.t b )=0.34 y L y 
(Q 2 ) = 1.99 GeV 2 \_y 



.n = Lj L 1 1 1 L ! 1 1 1 \1 : 1 1 N i — U i 1 1 1 

90 180 270 90 180 270 90 180 270 90 180 270 360 



(f> [dcg] 4> [dcg] 4> [dcg] <t> [dcg] 

Figure 5: The beam spin asymmetry Alu(0) measured by the CLAS collaboration [74] for 
a sample of 4 bins, all taken at an average —t value of 0.3 GeV 2 . For other notations refer 
to Fig. 3. 



3.3 Longitudinal target polarization 

This asymmetry is very similar to the beam spin asymmetry in that it is dominated by 
the DVCS-BH interference. The most important observable of this type is A^ nr ^ which is 
primarily sensitive to % as an inspection of Tab. 4 reveals. The comparison of the HERMES 
data [96] on this observable measured with a positron beam, with our predictions is made in 
Fig. 6 and, given the large experimental errors, reasonable agreement is observed. Only at 
small —t our result seems slightly larger in absolute value than the data by a fraction of the 
HERMES error bar. A surprisingly large sin 20 harmonic appears in the HERMES data, 
see Fig. 6. Theoretically this contribution should be heavily suppressed to the order we 
are working and indeed our prediction is very small. The reason for the large experimental 
value of Ay^ 2 ^ is unclear. For both A^} L parameters the figure also shows the contribution 
from H separately, confirming that this GPD represents a significant fraction of A^ n( ^. 

The CLAS collaboration has published data on A^ L [99] , integrated over a large kine- 
matical bin. In contrast to HERMES this collaboration uses an electron beam implying 
that the BH-DVCS interference contributes with opposite sign to Ay L than in the HER- 
MES data on A^ L - The sin0 and sin 20 harmonics are shown in Fig. 7 along with our 
predictions. In contrast to the HERMES data, the sin 20 harmonic is compatible with zero, 
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Figure 6: Top plots: The sin0 and the sin 20 harmonics of 



and the cos0 harmonics of A^ L . HERMES data are taken from Ref. [96] 



UL . Bottom plots: The cos 00 
the average 

kinematics is xb = 0.1 and Q 2 =2.46 GeV 2 . The dashed lines in the top plots represent the 
results with only the contributions of "H and the dotted lines in the bottom plots the BH 
contributions. Our results shown as solid lines with the shaded areas as the error bands, are 
evaluated at the specified kinematics corresponding to each bin in —t. For other notations 
refer to Fig. 3. 



however our prediction is small and positive, off by about 1.5 a. For the sin0 harmonic, 
our result and the CLAS data are in perfect agreement. 

HERMES has also measured the double-spin asymmetry for a longitudinally polarized 
beam and target. In contrast to the single spin asymmetries Ay L and A\ v the BH term 
contributes to A\ L . In fact it provides the dominant contribution to the cos 00 term 
and about half the magnitude of the cos0 harmonic, see Fig. 6. Thus, at the current 
experimental accuracy in these observables, no information about GPDs can be extracted. 
Our full results for A\ L are of course in agreement with HERMES data as can be seen 
from Fig. 6. 
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Figure 7: CLAS longitudinal target single spin asymmetry A UL [99] at the average kine- 
matics x B = 0.28, -t = 0.31 GeV 2 and Q 2 = 1.82 GeV 2 . The left plot shows the sin0 
moment versus —t whereas the right plot shows the sin 20 moment. In both plots, our 
results are shown as solid lines with the shaded areas as the error bands. 



3.4 Transversely polarized target 

The sin(0 — <f>g) harmonic of the transverse target spin asymmetry is especially interest- 
ing since it is generated by the squared DVCS contribution and it is mostly sensitive to 
Im('H£*), see Eq. (56) and Tab. 4. For given H it is uniquely sensitive to the CFF £. It 
therefore constitutes one of the few ways to constrain the GPD E. The valence part of 
E is rather well constrained by fits to the data of the nucleon form factors [36] and is in 
fair agreement with the data on Aut for p° electroproduction [52, 53]. As we discussed in 
Sect. 2.2 for p° production the contributions from E for sea quarks and gluons cancel to a 
large extent. In view of this the DVCS data on Aut are complementary to DVMP and es- 
pecially important because to LO QCD there is no contribution from E 9 and therefore the 
cancellation between the gluon and sea-quark contributions cannot happen. Hence, Aut 
for DVCS probes the sea-quark part of E for given GPDs H and -Eyai- We remark that 
^Ijt dvcs ls f° rce d to vanish for forward scattering by angular momentum conservation. 

The HERMES data [95] on the sin(0 — (/)$) harmonics of Aut is shown in Fig. 8. Also 
shown in Fig. 8 is A^^~^ s ^ cos ^. rpj^g observable receives separate contributions from ImH, 
ImS and IhiH (see Eq. (56) and Tab. 4) and even for E = H = reasonable agreement 
with experiment is achieved. In contrast to this observable ^ut bvcs would vanish if E 
was zero which is not the case experimentally. 

In order to check the normalization of the GPD E for a flavor symmetric sea we also 
display in Fig. 8 results for three different normalizations of E s : N s = and N s = ±0.155; 
for the latter value a positivity bound is saturated [55] (see Sect. 2.1.2). It can be seen 
from the figure that a negative E s seems to be favored. Considering the large errors of the 
data, a value of N s close to zero but negative is not excluded while a large positive value 
is apparently in conflict with experiment. 
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Figure 8: The transverse target spin asymmetries. The left plot shows ^uTDvefs whereas 
the right plot shows ^uti" * s ' )cos ^. Data are taken from Ref. [95], their averaged kinematics 
is %b = 0.09 and Q 2 = 2.5 GeV 2 . Our results are shown for the standard scenario (N s = 
—0.155 ) as well as for two other normalizations of E Bea (N s = and 0.155) and finally 
keeping only GPD H and setting the others to zero. For other notations refer to Fig. 3. 



3.5 Hall A cross sections 

The helicity-dependent ep — > epj cross section has been measured by the Hall A collab- 
oration at Jefferson Lab [73] at fixed ££=0.36. These data are extremely accurate and 
are therefore very demanding to theory. They have the advantage to cleanly separate the 
difference of the cross sections for opposite electron helicities, Act, from their sum, Ecr. The 
main contribution to the cross section difference is the sin harmonic fed by the BH-DVCS 
interference. From Tab. 4 we see that the most prominent contribution from the CFFs 
is IniH. The cross section sum, on the other hand, which corresponds to the unpolarized 
cross section, receives contributions from all three terms BH, DVCS and the BH-DVCS 
interference but the BH contribution is dominant. Its most prominent harmonics are cos 00 
and cos0. 

The comparison of the difference and sum of the Hall A helicity cross sections with our 
results at their highest Q 2 bin is shown in Fig. 9. The agreement is nearly perfect for Act 
which indicates that our parameterization of H seems to be adequate. This is in line with 
our findings for A^^ at HERMES kinematics, namely that our predictions agree well with 
the recoil data [98]. 

The cross section sum exhibits the expected dependence of a linear combination of the 
cos 00 and cos harmonics. There is however a discrepancy between theory and experiment 
of up to 30% in the strength, especially around = 180°, pointing to something missing 
in the real parts of the various CFFs entering this observable, in particular Re"H, see 
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Figure 9: Jefferson Lab Hall A helicity-dependent cross section data at and different t bins 
for xb = 0.36 and Q 2 = 2.3 GeV 2 . The top plots show the differences of cross sections for 
opposite electron helicities versus <j) whereas the bottom plots show the unpolarized cross 
section. Data are taken from [73]. The Bethe-Heitler contribution to the unpolarized cross 
section is represented by dashed lines whereas our full results are shown as solid lines with 
the errors as shadowed bands. 



Tab. 4. Note that by taking the ratio of the difference over the sum of cross-sections 
for Hall A, one recovers a result fully compatible with the CLAS beam spin asymmetry 
data in the same kinematical range. It is therefore no surprise that our approach also 
has difficulties reproducing the CLAS beam spin asymmetries since the issue apparently 
lies in the unpolarized cross section. The real part of H is also probed by A c ^ sn< ^ . Since 
KeH contributes to Xcr with a small coefficient as compared to the BH-term, a substantial 
increase of Re'H would be required in order to fit the Hall A Sa data. Whether such 
an increase is compatible with other DVCS and DVMP data remains to be seen. In any 
case, we would like to stress again that the GPD parametrization we are using has been 
tuned to much lower values of £. Thus, for instance, the addition of a D-term to the 
double-distribution might improve the agreement between theory and experiment. The 
exploration of this as well as other improvements of the GPDs is left for future work. 
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3.6 HI and ZEUS cross sections 



The HI [91, 92] and ZEUS [93, 94] collaborations published t-differential 7*p ->■ 7p cross 
section. The kinematics for these data is characterized by small xb, typically 1CT 3 — 10~ 4 , 
a photon virtuality that varies between 3 and 25 GeV 2 and large W, of order 100 GeV. In 
this kinematical range the GPDs at small skewness and small x control the DVCS cross 
section. The dominant contribution to the cross sections comes from the GPD H. The 
GPD E is suppressed, see Tab. 4. The contribution from H, although not suppressed as 
compared to H (see Tab. 4), is negligible here since H is much smaller than if; this is in 
particular the case for the sea quark contribution. We remark that for the leptoproduction 
of vector mesons, H does not contribute to leading-twist accuracy. Due to the large range 
of Q 2 in which the HERA data are available, evolution of the GPD H plays a decisive role. 
Using the scale-dependent parametrization of H introduced in Sect. 2.1.1 we evaluate the 
DVCS cross section at HERA kinematics to leading-twist accuracy and LO of perturbative 
QCD. As the factorization scale we choose hf — Q- Our results are compared to the HERA 
data in Fig. 10. Reasonable agreement with experiment is achieved within experimental 
errors and theoretical uncertainties. A similar observation has also been made in [31]. In 
this paper only the GPD H is taken into account and parametrized in terms of an SO (3) 
t-channel partial wave expansion. A combined fit of the GPD parameters to the HERA 
DVCS and DVMP data is performed. In contrast to the analysis carried through in [25, 26] 
(see Sect. 2.2) DVMP is also computed within the collinear factorization approach. This 
necessitates a GPD that differs from the one we are using in particular at the cross-over 
line £ = x, and that leads to very strong evolution effects. The quality of the combined fit 
to the HERA data on DVCS and DVMP performed in [31] is comparable to that of our 
results. 

4 Future Experiments 

The COMPASS collaboration will have a two-week DVCS run in 2012 with a 160 GeV 
muon beam, followed by a longer run in 2015 and later [101]. They plan to measure 
the t-slope of the ^-integrated photon electroproduction cross section, as well as specific 
charge and spin observables. The polarized muon beam is produced through pion decay, 
its polarization therefore changes sign when the beam charge is reversed, i.e. fi + and /i~ 
are polarized along opposite directions. COMPASS plans to measure mixed charge-spin 
(CS) cross sections differences and sums defined as follows: 

<S CS) u(0) = da* + da^ = 2da uu {l + A LUjI (0)) 
^cs,u(0) = da* -da^ = 2da uu (A hV ^ cs (4>) + A c (4>)) 

where the cross sections da h >*> e » are defined in Eq. (44). Predictions from our approach for 
these three observables are shown in Fig. 11 along with predictions made in Ref. [24]. 
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Figure 10: Differential DVCS cross section versus t for a set of Q 2 values and large W 
values ranging from 71 GeV at low Q 2 to 104 GeV at the highest Q 2 . Data are taken from 
Refs. [91, 92, 93, 94], where statistical and systematical errors are added in quadrature and 
normalization uncertainties were ignored. Our predictions are shown as solid lines with 
errors represented by shadowed bands. 
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Figure 11: Charge-spin asymmetry (left), cross section difference (middle) and sum (right) 
as a function of for one of the COMPASS-II kinematical bins at xb = 0.05, Q 2 = 2 GeV 2 
and —t = 0.2 GeV 2 . Our predictions are shown as solid lines with error bands, the 
predictions for the two scenarios given in Ref. [24] are represented by short dashed and 
dash-dotted lines. 

The CLAS12 collaboration plans to measure the beam and target spin asymmetries as 
well as cross sections for DVCS with the 11 GeV electron beam as soon as 2015 [102]. In 
Fig. 12 our predictions for the sin0 moments of the beam and target spin asymmetries 
^Lu m< ^ an d ^uL m< ^' evaluated with the help of the analogue of (54), are displayed for a 
typical kinematical bin accessible with CLAS12. The Hall A collaboration will also run a 
DVCS experiment very soon after the 12 GeV upgrade is complete [103]. Our predictions 
for the difference and sum of the helicity cross sections are shown in Fig. 12 for one of the 
planned kinematical settings. 

More detailed predictions can be obtained from the authors on request. 

5 Summary and outlook 

Since 2000 there has been a great wealth of measurements related to the DVMP and 
DVCS processes. Factorization theorems assert that these hard exclusive processes can 
be interpreted in terms of partonic degrees of freedom, provided the virtuality of the 
exchanged photon is large enough. In that case the partonic picture involves GPDs, which 
are universal, process-independent quantities. Up to now, both processes have mostly been 
studied independently of each other, and it is therefore an important check of consistency 
to test whether the GPDs used to describe one channel can be used to describe the other 
channel as well. Obviously universality is exploited at the most in a combined analysis of 
both DVMP and DVCS but this large-scale program is beyond the scope of the present 
article. 

A necessary and interesting first step in exploiting universality is the use of a set 
of GPDs extracted from an analysis of DVMP data from HI, ZEUS, E665, COMPASS 
and HERMES [26, 28, 55] in a then parameter-free evaluation of DVCS and the detailed 



29 




Figure 12: Upper left and right: the sin0 harmonics of the beam and target spin asym- 
metries A^j and ^ul^ versus — t for a typical bin accessible with CLAS12 [102]. Lower 
left and right: the difference and sum of the helicity cross sections versus in one of the 
kinematical bin of the Hall A DVCS experiment [103] . Our predictions are shown as solid 
lines with error bands. 
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comparison with the available ep — > ejry data. The various observables are computed within 
an approach in which the DVCS amplitudes are calculated to leading-twist accuracy and 
leading-order perturbative QCD while, similarly to Ref. [33], the leptonic tensor is taken 
into account without any approximation. We stress that DVCS is treated in the collinear 
factorization framework. This is consistent with the calculation of DVMP as described 
in Sect. 2.2 and hence, the GPDs extracted from DVMP can be used for the calculation 
of DVCS observables. In both processes the quarks are emitted and reabsorbed from the 
proton collinearly with the proton momenta. 

We observe very good overall agreement between our predictions and most of the HI, 
ZEUS and HERMES data but a less satisfactory description of the large xb (small W) data 
from Jefferson Lab, where power corrections might be large. We notice that the most recent 
HERMES measurements [98] using a recoil detector in order to achieve fully exclusive final 
states, give a significantly larger beam spin asymmetry (in absolute value). The recoil 
data are in perfect agreement with our results. Such dilution effect may be present in 
other HERMES observables as well and estimates of this effect would be highly valuable. 
The COMPASS, CLAS12 and Hall A collaborations have planned to take DVCS data in 
the coming years. In view of the good agreement with the available data demonstrated in 
our study we also give predictions for their main observables. 

Our study makes it clear that the set of GPDs we are using which is extracted from 
DVMP, describes DVCS data over a large kinematical range rather well in general although 
not perfectly in all details. This is partly a consequence of the different kinematic ranges 
of DVCS measurements and of the DVMP data used to constrain our GPD model. In 
addition, the GPD parametrization is of course an approximation. Thus, improvements of 
the GPD parametrization are required on the long run. In order to do so future data from 
COMPASS and JLabl2 will be of help. Possible improvements may include the use of more 
recent versions of the PDFs, the proper scale dependence of the GPDs, an updated form 
factor analysis [54] , eventual modifications of the profile functions of the double distribution 
parameterization and allowance for a non-zero D-term. 
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